Inelastic resonant tunneling through single molecules and quantum dots: 
spectrum modification due to nonequilibrium effects 
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Resonant electron transport through a mesoscopic region (quantum dot or single molecule) with 
electron-phonon interaction is considered at finite voltage. In this case the standard Landauer- 
Biittiker approach cannot be applied. Using the nonequilibrium Green function method we show 
that due to a nonequilibrium distribution function of electrons in the mesoscopic region, the inelastic 
scattering rate and spectral function of the dot become functions of the voltage and have to be 
calculated self-consistently. 



During the past several years nonequilibrium quantum 
transport in nanostructures and, in particular, transport 
through mesoscopic tunneling structures, quantum dots 
(QD) and single molecules, is in the focus of both experi- 
mental and theoretical investigations because of possible 
electronic device applications. In this paper we recon- 
sider one of the basic problems, namely resonant tunnel- 
ing through a mesoscopic region with electron-phonon 
interaction, placed between macroscopic leads. We want 
to show that the nonequilibrium distribution function 
of electrons in the dot changes significantly the spectral 
function of the dot. 

In the absence of interactions and in the limit of small 
voltage the current through a QD is described by the 
Landaucr-Buttikcr formula 



J = 



2irh 



T(e)[f(e-e^ L )-f(e-e VR )]de, (1) 



where T(e) is the transmission coefficient, (fiL(R) is the 
electrical potentials of the left (right) lead. In the case 
of one resonant level with energy eo, the transmission 
coefficient is given by the Brcit-Wigner formula T(e) 



r 2 / [(e 



eo 



r 2 ] , where T is the level width produced 



by coupling to the leads. The corresponding spectral 
function of the dot, current-voltage curve and differential 
conductance are shown by the thin lines in Figs.^ HI 

At finite voltage the current can be described by |J1J 
with a voltage-independent transmission coefficient only 
in some special cases, e.g. for energy-independent cou- 
pling to the leads (wide-band limit). In the presence 
of interactions (e.g. internal Coulomb interaction and 
electron-phonon interaction ) T(e) becomes temperature 
and voltage dependent. 

The problem of inelastic resonant tunneling of elec- 
trons coupled to phonons was first considered in Rcfs. 
Ill Li Li Phonon satellites are formed in the spectral 
function and in the differential conductance. Typical re- 
sults are shown by dashed lines in Fig. ^ and Fig. [21 Re- 



cently this problem attracted attention after experiments 
on inelastic electron transport through single molecules 
Li Li lil Li III B- New theoretical treatments were pre- 
sented in Refs. [ljj, [TJ Hj, U&, LLil In our paper we study 
the influence of the nonequilibrium electronic distribu- 
tion function on the phonon self-energy. We find a sig- 
nificant shift and broadening of the phonon-satellites in 
the electronic spectral function. To our knowledge this 
effect has not been considered yet. 

We use the Nonequilibrium Green Function (NGF) 
method |la, |l6|, which now is a standard approach in 
mesoscopic physics and molecular electronics [ijj [lj, UM 
|2fJ . The advantage of the NGF formalism is that it can 
be successfully applied to a variety of systems and prob- 
lems, that it is initially exact, and many pow erful ap- 
proximations can be derived from it (see | 2(1 |21| and 
references therein). 

We describe the electronic states of the mesoscopic re- 
gion (QD) by a set a single-electron states \a) with ener- 
gies e a . These are coupled to the free condcuction elec- 
trons in the leads by the usual tunneling Hamiltonian. 
Furthermore the dot-electrons are coupled to vibrational 
modes. We do not consider Coulomb interaction in the 
dot, to avoid further complications, such as Coloumb 
blockade and Kondo effect. The system is then described 
by the following model Hamiltonian: 

H = H D + H L + H R + H T + H ph , (2) 

with the dot-Hamiltonian 

H D = Y,(t a + e(p D (t))did a , (3) 

a 

the Hamiltonian of the right (i = R) and left (i = L) 

lead 



H l=L ,R = ^(ejfecr + e(pi{t))c] ka .c 



ika i 



ho 
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and the tunneling Hamiltonian 



(4) 



(5) 



Finally H p h describes phonons and clcctron-phonon cou- 
pling 

H p e t = Y, hbJ A a i + J2J2 V(a- 9 + a\)did a , (6) 



where fioj q is the phonon energy. The quantities ipi (t) 
and <fD(t) are the electrical potentials of the leads and 
the dot The latter will be set to zero in the following. 

The current in direction from the left (i = L) or right 
(i = R) contact to the dot is found to be 



, ^ 2e 
Ji(t) = -r- Re 



Vikcr,aG alka (t,t) 



(7) 



where we define the lead-dot lesser Green function 
Ga.tko-fa^a) = i(4fc<7(*2)da(*i)), while the charge of 
the QD is defined by the dot lesser function G5g(ii, fe) = 

i(dl(t 2 )d a (h)y 

The equation of motion for the lesser Green func- 
tion G < contains coupling to the retarded and advanced 
Green functions G Rt . In order to obtain a simple Dyson 
equation it is convenient to combine these functions to a 
matrix in Keldysh space [id, I20L |21| 



G A 



(8) 



where each component again is a matrix G,,.,/^!,^) 
within the states of the system including all the states 
of the leads (77 = ika) and the states of the dot (r/ = a). 
For this matrix Green function the noncquilibrium 
Dyson equation can be written in the usual way, in dif- 
ferential form 



'sr-*< fa > 



G 



{to) 



S, 



(9) 



where 6 = IS^'S^i — t 2 ), an d / is the unit matrix in 
Keldysh space. H(t) = (H L + H R + H D + H T )I is 
the single-particle Hamiltonian. The curly brackets de- 
note a matrix multiplication and convolution in time: 

{AB} m ,(t 1 ,t 2 ) =-- J2 I dt 3 A ril {h,t 3 )B w (t 3 ,t2). The 



E 

7 



self-energy S describes interactions in the dot. The re- 
tarded and advanced Green functions and the spectral 
function A = i (G R — G A ) describe quasiparticlc exci- 
tations. The lesser function G < contains in addition a 
quantum distribution function. 

In the absence of interactions in the leads (besides the 
tunneling) one can derive the following exact expression 
for the lead-dot function [l7ll20|: 



Gq,. 



ika 







V$,iko J GapGika f - , (10) 



where G a p is the full dot Green function, while Gika is 
the bare Green function of the leads, which is diagonal 
in the quantum numbers ika. 

Substituting i|lU|) in the general equation © we obtain 
the following basic equation for the dot Green function 
G(ii, t 2 ) = G a /3(ti,t 2 )- 



d-H 



dti 



G-{±G}=I6(h-t 2 ), (11) 



withH = H al3 = e a S a/j . HcreS(e) = sf ) (e)+S^ T) (e) + 
S^ p/l ^(e) is the total self-energy of the dot composed of 
the tunneling self-energy 

S j=i,iJ = S ja/3 = 22 \ V jka,aGjkaVjka,0j ■ (12) 

kc 

and the self-energy Xw'' 1 ) = fy£„ ' describing interactions 
with phonons. 

In the following we consider stationary trans- 
port produced by a time-independent voltage. In 
that case the Green functions depend only on 
the time difference. After Fourier transformation 
(G(e) = / G(£i - t 2 )e"(* 1 -* 2 )d(ii - t 2 )) one obtains 



(eI-H)G(e)-£(e)G(e)=I, 



(13) 



m 



For the current we obtain the well-known expression 



Ji=|/^Tr{r i (e-e^)(G<(e)- 



(14) 



+/0(e-e^)[G*(e)-G^(e)])}, 
with the level-width function 
r i= i(ij)(e) =r M/3 (e) =2ir^2Vika,0V^. (rta 5(e-tika)- 

her 

Equations for the retarded and advanced functions fol- 
low from the diagonal part of (|13fl 



(e e a )G a p 2_^ ^q 7 G 7/ 3 — 5 a f: 



(15) 



and the equation for the lesser function (quantum kinetic 
equation) follows from the off-diagonal part of (|13|l com- 
bined with its conjugate G(e) (el — H) — G(e)S(e) = I: 



(ep - e a)G Q/3 - 2_^ ( ^a 7 G 7/3 + S Q7 G 



'7/3" 



-G Q7 S 7/3 - G Q7 S 7/3 ) — 0. 



(16) 



In the following we will concentrate on the tunneling 
through one level a only. Then l|16|l simplifies further: 

(E£ - S^)G< - (G R G A )X> = 0. (17) 



For the evaluation of the different terms in this equa- 
tion we make the following ansatz for the lesser Green 
function of the dot 



G<{e) 
A(e) = : 



= W(e)/(e). 



(18) 
(19) 



introducing a general nonequilibrium distribution func- 
tion /(e). For the Green functions of the leads in equi- 
librium we may use the Fermi function /(e). 

For the retarded tunneling self-energy T, R ^ TS) one ob- 
tains 



,fl(T) 



d 3 fc 



\V ik J 



(2tt) 3 



e fc - ey»i + iO 



A(e-eipi)--Ti(e-eifi), 

(20) 




where A is the real part of the self-energy, which usually 
can be included in the level energy e Q , and T describes 
level broadening due to coupling to the leads. For the 
corresponding lesser function one finds 



,<(T) 



iTi(e-&fi)f (e - eipi). 



(21) 



In the following we will neglect the energy dependence 
of the level- width- function T. Then from (|17J) we obtain 
the following kinetic equation for the nonequilibrium dis- 
tribution function /(e): 

Tl [/(e) - /°(e - &p L )] + T R [/(e) - /°(e - &p R )] + 



iE<Cpfc)( e ) + i h: R (P h )( e ) - S A (P' l )( e )) /( e ) 



0. 



(22) 

Finally we have to evaluate the self-energy T,^ ph ' due 
to the interaction with phonons. In order to demonstrate 
nonequilibrium effects most clearly we consider only one 
phonon mode with frequency ujo. In the standard self- 
consistent Born approximation, usin g K eldysh technique, 
one obtains for the self-energies (2f 



■£R(ph) 



2 
2 



J^(G R _ e ,D« + G?_ e ,D*), (23) 



'111. 

2vr 



XAiph) = ™_ J ^_ ( G A_ e/D K + G K_ eiD A^ ^ (24) 



de' 



Z<(ph)=i\* I —G< f ,D<. 
2tt e_e 



(25) 



where G K = 2G< +G R -G A is the Keldysh Green func- 
tion. 

In the following we do not consider phonon renormal- 
ization and polaron effects, and use bare phonon Green 
functions: 



£ R (e) = 



1 



1 



(26) 



e — ljo + iS e + ljq + i5 ' 
D<(e) = -2m [(N + l)5(e + lu q ) + N Q 5(e - w )] (27) 



Figure 1: (Color online) Spectral function of the dot with en- 
ergy level eo = 2. Thin line: no interactions, dashed line: in- 
teraction with phonon, electronic distribution in equilibrium, 
thick line: with phonon, including nonequilibrium effects at 
high voltage 



> 




Figure 2: (Color online) Current-voltage curves and differen- 
tial conductance (eo = 2). For the legend of the curves see 
Fig. 1 



where iVo = l/(e w °/ T — 1) is the equilibrium phonon 
population, (nonequilibrium effects in the phonon dis- 
tribution was considered recently in |l4l]). However, 
we include nonequilibrium effects in the electron dis- 
tribution function while calculating the self-energies. 




Figure 3: (Color online) Differential conductance at small 
voltages for different positions of the dot energy level eo = 
(circles), eo — 1 (squares), and eo = 2 (triangles). 



These self-energies enter the spectral function A(e) = 
i(G^(e) — G„ (e)) of the quantum dot, where G^(e) can 
be found from (fT5|l 



G£(e) = 



1 



,R(ph) 



+ i(T L +T R )/2 



(28) 



The real part of the self-energy leads to an energy shift 
in the spectral function, which depends on the applied 
voltage. To take this effect into account properly, we per- 
form self-consistent calculations of the distribution func- 
tion /(e) ill'L'li and self-energies (J23J1 - H25J) a t finite voltage. 
The spectral function is then used to calculate the cur- 
rent voltage curves. Combining Jl and J R (see [lT]]) the 
expression for the current can be written as current 



J 



rv.r 



L*-R 



2-kTl T r + Ti 



deA(e)[f°(e-e<p L )-f°(e-e<p R )}. 

(29) 

It looks as simple as the Landauer-Biittiker formula Q), 
but it is not trivial! The spectral density A(e) now de- 
pends on the distribution function /(e) and hence the 
applied voltage, ip^ = —fR = V/2. Results for the 
current-voltage curves are presented in Fig.[3] by thick 
lines. 

If the spectral function A(e) is calculated with help of 
the equilibrium distribution function /°(e), which does 
not depend on the voltage, the usual phonon satellites in 
the spectral function and differential conductance are ob- 
tained as shown by the dashed lines in Fig.^and Fig-El 
But if one takes into account the noncquilibrium distri- 
bution function /(e), the spectral function changes with 
increasing voltage from the equilibrium one to the func- 
tion shown as thick line in Fig.^for a voltage V > eo. 
Correspondingly, the phonon resonances in the differen- 
tial conductance are shifted and suppressed (Fig.0 thick 
line). 

Finally in Fig.[3]we show the differential conductance 
for different positions of the dot level eg with respect to 



the Fermi energy of the leads (defined in the absence of 
an applied voltage). For small dot energies the phonon 
peak appears at zero voltage and broadening is weak. 
This property can be used for experimental investigation 
of inelastic effects by tuning the dot level with help of a 
gate voltage. 

In conclusion, we investigated inelastic resonant trans- 
port through a mesoscopic region. As an example, we 
considered resonant electron tunneling through a single 
level coupled to a single phonon mode. This model, al- 
though quite simple, is popular in connection with trans- 
port through vibrating molecules. The self-consistent 
treatment of the electronic nonequilibrium distribution 
function shows that the main peak and phonon side- 
bands in the differential conductance are shifted consid- 
erably and are essentially broadened as compared with 
their equilibrium position. 

We thank K. Richter and G. Cuniberti for valuable 
discussions. 
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